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SOBOLEV-TRACE INEQUALITIES OF ORDER FOUR 


ANTONIO G. ACHE AND SUN-YUNG ALICE CHANG 


Abstract. We establish sharp trace Sobolev inequalities of order four on Euclidean 
d-balls for d > 4. When d = 4, our inequality generalizes the classical second order 
Lebedev-Milin inequality on Euclidean 2-balls. Our method relies on the use of 
scattering theory on hyperbolic d-balls. As an application, we characterize extremals 
of the main term in the log-determinant formula corresponding to the conformal 
Laplacian coupled with the boundary Robin operator on Euclidean 4-balls. 


1. Introduction 


In this paper we derive some sharp Sobolev trace inequalities on Euclidean d-balls 
of order 4 with the {d — l)-sphere as boundary. Sobolev inequalities for critical 
exponents play an important role in problems in conformal geometry. One example 
is the sharp Sobolev inequality of order 2 on the standard sphere {S"',gsn) for n> 2 
which takes the form 


( 1 . 1 ) 



g-2 

n ■ vol(S'”) 


/ |V/pda 


1 

vol(S'"') 



I/pda, 


where da is the Lebesgue measure induced by (75^, vol(S'”) is the volume of S'" with 
respect to V is the gradient operator on {S^,gs^) and q = Inequality (II.Ih 
has been crucial in the study of the Yamabe problem on closed manifolds. We remark 
that when n = 2 , fll.ip takes the form of Moser-Onfori inequality 


which is fundamental in the problem of prescribing Gaussian curvature on the sphere 
(see f |Mos71] . |Ono82] . see also |OPS88a( IOPS88b( IOPS88c] )b On manifolds with 
boundary, a class of inequalities analogous to fll.ip and fll.2p is the class of Sobolev 
trace inequalties. A classical example of these inequalities is the following: consider 
the unit ball B'^ on Euclidean space with the sphere as boundary and let / be 
a function in If n is a smooth extension of / to we have the inequality 

/ l/l^daV”' < / |Vnpdz + ^^/ fda 

J S‘^-^ J J B‘i 2 Jsd-I 

when d > 2. 


(1.3) lL2voi(s-'-')jii 
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When d = 2, inequality fll.3p becomes the classical Lebedev-Milin inequality |LM51j 

(1.4) log ( — (f e'^da] <— [ \'Vv\'^dx — (f fda, 

VTT Jsl J 4:71 Jjy 71 Jgl 


where now Vn is the gradient of v with respect to the Euclidean metric on the disk 
= Inequality fll.Sp has been derived by Escobar [Esc92] 

and applied to study the Yamabe problem on manifolds with boundary; while the 
Lebedev-Milin inequality fll.dp has been applied to a wide variety of problems in clas¬ 
sical analysis, including the Bieberbach conjecture |dB85] 1 and by Osgood-Phillips- 
Sarnak f |OPS88^ IOPS88bl IOPS88cj i in the study of the compactness of isospectral 
planar domains. 

In this note, we consider extensions of some of the above results to order 4. The 
order here refers to the order of operator involved in the derivation of the inequalities. 
For the inequalities cited above, the operators which describe the Euler-Lagrange 
equation of the extremal functions involved in the proof of the inequality are either 
Laplace operator or the conformal Laplace operator, both of order 2. In recent years, 
the role played by these second order operators has been vastly extended to operators 
of order 4 (e.g. the bi-Laplace operator or the Paneitz operator of order 4, which we 
will briefly discuss in section [2]). For example, the Sobolev inequality which prescribes 
the embedding of to on R” for n > 4 has a fourth order analogue on S^, 
namely 


(1.5) 


a 


vol(S^ 


I/I--4 da 


< 




fPifda 


where Cn = 2)(n+2){n 4 ) ^ order Paneitz operator on dehned 

as 


Pa = 



n 



n 


■) 


(see |Bra951 Theorem 1.1]) where A is used to denote the Laplacian with respect to 
the metric (see |Bec931 Theorem 6]. When n = 4, the inequality becomes an 
inequality of Moser-Onofri type just as in fll.2p . and takes the form 

which is a special case of Beckner’s log-Sobolev inequality, see |Bec93] , also see Chang 
f [Cha96j . Lecture 5). As mentioned above, the main results in this paper are ex¬ 
tensions of inequalities fll.3p when d > 2 and fll.4p above to sharp Sobolev trace 
inequalities of order 4. We will adopt the following notational convention 


Convention 1.1. We will sometimes use go to denote the Euclidean metric on R'^. 
The Euclidean metric is also commonly denoted as |da:p. 


We start by stating our result for dimension d > 4 
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Theorem A. Given f G with d > 4, suppose v is a smooth extension of f 

to the unit ball which also satisfies the Neumann boundary condition 


(1.7) 


d , 

on 


(^- 4 ) 

4 


Then we have the inequality 






d-4 

d-1 


( 1 . 8 ) 


< [ (^g^vfdx + 2 if |V/|^<itr + bdf Ifl^da, 

J B'i Jsd--^ J s<^-^ 


where Cd = = d{d^ andV f is the gradient of f with respect to the round 

metric gsd-i- Moreover, equality holds if and only if v is a bi-harmonic extension of 
a function of the form fz^H) = c|l — (zq, where c is a constant, f G zq is 

some point in and v satisfies the boundary condition fll.7p . When / = 1, equality 
in 01.81) is attained by the function v(x) = 1 + 4^(1 — \x\‘^). 


The corresponding theorem for dimension 4 is the following inequality which is a 
generalization of the classical Lebedev-Milin inequality 01.4^ : 

Theorem B. Given (p G C°^{S^), suppose w is a smooth extension of (p to the unit 
ball B'^. If w satisfies the Neumann boundary condition 

(1-9) = 0, 

then we have the inequality 

( 1 . 10 ) + 

Moreover, equality holds if and only if w is a bi-harmonic extension of a function of 
the form (fzoiO = ~ log |1- ~ (^Oi 01 + 0 where c is a constant, f G S^, Zq is some point 
in B^ and v satisfies the boundary condition 01.91) . 


We remark that the method that we used to discover and establish the above 
sharp inequalities is non-traditional in the sense that we hrst derived such inequalities 
for some metric g* on B^ which is in the conformal class of the Euclidean metric 
go = and then we derived the desired inequalities in Theorems lAl and iBl through 
conformal covariant properties between the 4th order “Paneitz operator” on g* and 
the bi-Laplace operator on the flat metric go = \dx\^. The choice of the metric g* 
in turn comes from the consideration of the hyperbolic metric 5 ^+ on B'^ and the 
connection between g* and g^ is established through the scattering theory. In section 
2 below we will explain this connection and why the choice of g* is a “natural” one for 
the Sobolev inequalities in Theorem and [B] above. We remark that when we apply 
the same procedure to the second order conformal Laplacian on B'^, the g* metric we 
obtain agrees with go. In this sense, our choice of the metric g* is natural. 
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We also point out that once the sharp inqualities in the statements of Theorem A 
and B are formulated, it is not surprising that one can find more direct proofs of these 
two theorems, but we would like to record here our original approach to discover such 
sharp inequalities. We also believe that our method may lead to the formulation of 
other sharp inequalities. 

1.1. Application to Polyakov-Alvarez type formulas. We conclude by provid¬ 
ing an application to the study of the extremals of the functional determinants for 
conformally covariant operators on manifolds of dimension four with boundary. 

First we recall the classical Polyakov-Alvarez formula for surfaces with boundary: 
if (M, dM, g) is a surface with boundary and we let denote the metric 

gw = e^^g, 

where w is a. smooth function on M, then there is a formula relating the functional 
determinants associated to the Laplace Beltrami operators Ag^ and Ag and this 
formula is written entirely in terms of integrals of local invariants of (M, g) and of 
dM. More precisely, we have 

( 1 . 11 ) - ^ kgUds)g^ - kg{ds)g^ . 

Here Kg represents the Gauss curvature oi g, kg is the geodesic curvature of dM with 
respect to g and {ds)g is the line element on dM with respect to g. In the celebrated 
work of Osgood-Phillips-Sarnak |OPS88a[IOPS88b[IOPS88cj it is shown that there is a 
relationship between the Lebedev-Milin inequality fll.dp and the study of extremals of 
the ratio w — ?■ F[w]. More precisely, if in fll.lip we isolate the functional determinant 
coefficient 

(1.12) Fi[w] = [ \Vw\ldVg + [ wKgdVg+(f wkgdsg\ 

then Fi [tc] is conformally invariant and it is proved in |OPS88b] that inequality fll.4p 
can be applied to show that Fi [tc] has a definite sign when the background metric is 
the Euclidean metric. 

A fundamental point in the derivation of (II.lip is that the Laplace Beltrami oper¬ 
ator Ag has a conformal invariance property in dimension 2, however, this conformal 
property fails in higher dimensions. A result extending fll.lip for functional deter¬ 
minants of conformally covariant operators on compact manifolds with boundary in 
dimension 4 was obtained by Chang and Qing in |CQ97a| based on an earlier prelim¬ 
inary version of the formula by Branson and Gilkey |BG94] . In |CQ97a| , Chang and 
Qing considered the Conformal Laplacian operator Lg given explicitly in dimension 
4 by 
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and they coupled Lg with a boundary operator, namely the Robin operator dehned 
as 


an 


-H, 


where n is the outward normal to the boundary dM and H is the mean curvature of 
dM with respect to g. The pair {Lg, iBg) and its conformal properties are fundamental 
in the study of the Yamabe problem with boundary (see for example |Esc92] b Let 
{M'^,dM,g) be a compact manifold with boundary and let g^, = e^'^g where tc is a 
smooth function dehned on M. For the pair, {Lg,^g), it is proved in |CQ97a| that 
there is a Polyakov-Alvarez formula of the form 


(1.13) 


/ det(L,..S,.) \ 

t det(L„»,) ) 


11 

'^-f,li[g,w], 

i=l 


where the terms li[g, w] for z = 1 ,..., 11 are computed in terms of integrals of local 
invariants of M and dM with respect to the metrics g and g^- Moreover the second 
term l 2 [g,w] is conformally invariant and plays the same role as Fi[t(;] in fll.lip 
(where is given by fll.l 2 p i. We will show an explicit expression for l 2 [g,w] in 

Section 6 (see equation fIb.Sp L This is why we refer to \ 2 [g-,w\ as the main term in 
the functional determinant expansion fll.lSp . In ( |CQ97a| , Corollary 3.4), under two 
additional geometric assumptions, one can show that 12 ( 5 ', tc] has a dehnite sign on 
the model case {B^,S^,go). 

As an application of the sharp inequality in Theorem [HI we will here show that the 
metric g* is a minimizer of the functional I 2 over the class of functions with prescribed 
boundary value and Neumann boundary condition. More precisely, letting Cg, and 
Cg, be the following classes of functions 

(1.14) Cg,= !^we C^{B^) : w\ss = if, -^w\ss = o| , 
and 

(1.15) (?,. = L e ^“(B*) : Xis» = V, = -1 j . 


Observe that any function y in Cg, can be written as y = w + p where to G C<^. We 
will prove that the main term in the expansion fll.ldp is non-negative for the metric 
g* and we also describe in detail the equality case: 

Theorem C. Given Lp G C°°{S^) with e^'^da = lA^I, then for all w G Cg,, we have 
the following inequalities 

( 1 ) l 2 [p*,tc] > 0 with equality if and only if p>{i) = — log|l — ^ where 

Zq G B"^ is fixed, f ^ S^, c is a constant and w is a bi-harmonic extension of ip 
in Cg,. It follows that g* is a global minimizer of I 2 for functions in the class 
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(2) l 2 [ 5 'o, x] ^ 0 where x ^ Equality is attained for (f = — log |1 — {zq, 01 + c 
where zq G is fixed, f G S^, c is a constant and x 'Is a bi-harmonic extension 
of ip in C^. 


Remark 1.2. In |CQ97a CQ97b| an inequality similar to the one stated in Theorem 
iBlwas proved under very restrictive conditions. See for example |CQ97b Lemma 3.4] 


In a forthcoming joint work of the authors together with Jeffrey Case, we will give 
some further applications of the lower bound in Theorem ICl and also extend the study 
to extremals for the log-determinant functional (ll.lJh . 


1.2. Acknowledgements. The authors would like to thank Jeffrey Case for many 
helpful discussions. 


2. Adapted Metrics 

As we have mentioned in the introduction, one of the key steps used in the proof of 
the inequalities in this paper is the consideration of a special metric called “adapted 
metric” which we will denote by g*, which lies in the conformal class of the flat 
metric Qq = Idxp on the Euclidean unit ball B'^ and is with totally geodesic boundary. 
The consideration of the metric g*, introduced in the work of Case-Chang |CC14] . is 
natural from conformal geometry point of view, and via the connnection of the metric 
as a preferred conformal compactihcation (as we will explain below) of the hyperbolic 
metric on B'^. 

We now recall the dehnition of the g* metric in a general setting and will derive its 
specihc form in the special setting of B'^ with the sphere as boundary, and when the 
corresponding boundary operator is of order three. 

2.1. GJMS operators, scattering theory. We hrst recall the dehnition of the 
GJMS operators via scattering theory |GZ03] . A triple M", gjf) is a Poincare- 

Einstein manifold if 

( 1 ) is (diffeomorphic to) the interior of a compact manifold with boundary 
dX = M”, 

( 2 ) 5 f+) is complete with Ric( 5 f+) = —ng^, and 

(3) there exists a nonnegative r G C°°{X) such that r“^(0) = M”, dr 7 ^ 0 along 
M, and the metric g := r^gf^ extends to a smooth metric on X"'~^^. 

A function r satisfying these properties is called a defining function, since r is only 
determined up to multiplication by a positive smooth function on A, it is clear that 
only the conformal class [h] := [s'Itm] on M is well-dehned for a PoincarWEinstein 
manifold. 

Given a PoincarWEinstein manifold M"’, 5 f+) and a representative h of the 

conformal boundary, there exists a unique dehning function r, called the geodesic 
defining function, such that, locally near M, the metric g^ takes the form g^ = 
(dr^ -|- hr) for hr a one-parameter family of metrics on M with h^ = h and having 
an asymptotic expansion involving only even powers of r, at least up to order n (cf. 

[FGT^IGIMIIL^^ L 
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It is well-known (see |GZ031 IMM87j for more general statements) that given / G 
C°°{M) and s G G such that Re(s) > |, s ^ f+hIo 5 where Nq = NU{0}, and s{n — s) 
is not in the pure-point spectrum crpp(—of —the Poisson equation 

( 2 . 1 ) —lS.g^u — s{n — s)u = 0 in A, 
has a unique solution of the form 

( 2 . 2 ) u = Fr'^-^ + HF for F,H (A) and F\m = /• 


Indeed, F has an asymptotic expansion 

(2-3) F = /(o) -I- /( 2 ) 7 ’^ -I- /(4)7’^ + • • • 

for /(o) = / and all the functions f{ 2 tj are determined by /. The Poisson operator V{s) 
is the operator which maps / to the solution u = V{s)f, and this operator is analytic 
for s(n — s) ^ (Tpp(—Agp). The scattering operator is dehned by S{s)f = H\m- This 
dehnes a meromorphic family of pseudodifferential operators in Re(s) > Given 
7 G (0, ^), Graham and Zworski ( |GZ03] i dehned the fractional GJMS operator P 27 
as the operator 

(2.4) P2^f := d,S + t ) / for d-, = 2^’ 

For 7 G No, this dehnition recovers the GJMS operators |GJMS92] . For 7 G (O, |), 
Graham and Zworski showed that P 27 is a formally self-adjoint pseudodifferential 
operator with principle symbol equal to the principle symbol of (—A)^, and moreover, 
if h = is another choice of conformal representative of the conformal boundary, 
then 

(2.5) P 2 ,/ = e-—"P 2 , (e—V) , 

for all / G C°°{M). Together these properties justify the terminology “fractional 
GJMS operator.” 

We adopt the convention that for 7 G (O, |), the fractional Q-curvature Q 27 is the 
scalar 

(2.6) Q,-, := -^F„(l). 

In particular, we emphasize that this dehnition produces a well-dehned invariant in 
the critical case 27 = n, with the corresponding Q curvature a generalization of Bran¬ 
son’s Q-curvature (in the case when 27 = n = 4. see [Bra93a] . |GZ03] i. 


In this paper we will only deal with the special case when 7 = | and the correspond¬ 
ing GJMS operator P 3 of order 3. We will also treat P 3 as the “boundary operator” 
of the correponding GJMS operator. In the case 27 = 4, the GJMS operator P 4 is 
also called Paneitz operator and is an elliptic operator of order four dehned on the 
manifold X'^ (see |Pan08] I given by 


(2.7) 


(P 4 ), = (-A,)2 + ((4A, - (d - 2)Jgg)iX; •)) + 


d-4 


(Q 




2 
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where 6 denotes divergence, V denotes gradient on functions, Ag the Schouten tensor 
A = ^^(Ric — Jgg), Jg = 2 {d-i) ^ 9 ' ^9 scalar curvature of the metric g and 

is the Q-curvature 

( 2 - 8 ) {Q4)g = -^gJg + “ ‘^\^g?g- 

We remark that in general, properties and formulas for P 2 k and Q 2 k have been re¬ 
cursively derived 1 |Juh09] . |FG13] h The conformal invariance property of (P 4 ,Q 4 ), 
which is a special case of fl2.5p is best expressed as 

d+4 

(2.9) 

for all smooth function U dehned on manifolds of dimension d > 5, where g E [g] 

4 

denotes the metric g = {'tp)'^-'^g. and 

( 2 . 10 ) (Pt)iU = e-‘^Pt),(U), 

for all smooth function U dehned on manifolds of dimension d = 4, where g E [g] 
denotes the metric g = e^'^g. 

2.2. Adapted metrics when 7 = |. In the earlier paper of Case-Chang [CC14j . 
we have introduced a special metric called “adapted metric” to compactify a given 
conformally compact Poincare Einstein manifold 5A, 5 f+). The metric was 

introduced for any parameter 5 = ^ -|- 7 for 7 E (0, |) and for s = n when n is an 
odd integer. For such an s, we consider the Poisson equation fl2.ip with Dirichlet 
data / = 1, and denote the solution of the equation by Vg- By a result of J. Lee 
|Lee95j . if one assumes further that the Yamabe constant of the boundary metric h 
is positive (this happens when the scalar curvature of some metric in [h] is positive), 
then there is no point specturm of — and hence Vg > 0 on X. Thus we can 

take yg := (vg)^^ as a dehning function, and gg = ylg+ is dehned as the “adapted 
metric”. In the limiting case, when s = n and n is an odd integer, i.e. when 7 = f, 
the adapted metric as appeared in the literature in the work of Feherman-Graham 
f |FG02j . see also |YKG06] 1 and is dehned as follows: 

2 X d , 

( 2 . 11 ) r =- — \g=nVg, 

ds 

then r satishes 

( 2 . 12 ) -Ag+r = n, 
and the adapted metric is dehned as g* = e^'^g+. 

The adapted metric for all s satishes many good curvature properties stated in 
terms of “smooth metric with measures”. Here we will just cite the properties when 
s = ^ -|- I satished by the metric g* := gn_^_3 which we will use in the rest of this 
paper. 

Proposition 2.1 (Properties of g*). See Lemma 6.2 and Lemma 7.6 in [GG14] . 

(i) It is a metric with totally geodesic boundary. 
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(ii) Qii.9*) = 0 . 

(iii) We define the energy of U with respect to the Paneitz operator (^ 4)3 as the 
integral guantity obtained by dropping the boundary terms when integrating by 
parts of the integral fj^a{{P 4 )g*U}UdVg*, i.e., 

(2.13) Efig*)[U]= [ {AgMf - {4Ag* - {d-2)Jg.g*){VgM,Vg*U)dVg. 

Jxd 

For all smooth function f defined on dX^, we have the identity 

(2.14) W(g-)\U,] = / PUfda,. - ^ / Qz{a’)fda,.. 

^ Jax ^ JdX 

where Uf denotes the (unigue) solution of the eguation: 

(2.15) {Pi)g.Uf = Q on 

(2.16) Uf = f on dX, 

f)TT 

(2.17) —^ = 0 on dX. 

OUg* 

2.3. Explicit formula for g* metric on , gn)- In this section we derive 

explicit formulas for the metric g* in the model case of gn), which happens 

to be computable. 

Proposition 2.2. On the model case {B'^, gn), he., for g^ = p~‘^\dx\‘^, p = 
we have 

4 

(i) When d>5,g* = {i/j)'^\dx\‘^, where 

fi = l + P- 

(ii) When d = A, g* = e‘^^\dx\'^ = e^^~^^^^'>\dx\'^. 

Remark 2.3. The metric g* in (ii) is called Fefferman-Graham metric and has been 
described in several previous papers including |FG02l lYKCOGi ICC14] . The Fefferman 
Graham metric is constructed via scattering theory and satisfies several important 
properties, including 

• is totally geodesic under g* 

• The Q-curvature of g* is zero 

• Rg* > 0. 

2 

Proof, (i) Denote d = n + l, and s = | + |, by our definition, g* = Vs~‘’gn, where gn 
is the hyperbolic metric on B^ and Vg is the unique solution of the Poisson equation 
with Dirichlet data / = 1. 

(2.18) — Ag^Vs — s{n — s)vs = 0 in R 

We now recall gu = P~‘^go, where p = ^ and go = \dx\'^ the fiat metric on B^. 

Using this notation, we find that g* = fi'^-'^go and where = VsP ^ = VsP^~^. One 
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of the main observations in the paper of Chang-Gonzalez f [CGllj ) is that in the 
compactihed metric go = p^gn, equation fl2.18p is equivalent to the equation 

(2.19) + = 

where a = 1 — 27 = —2 in our case. The error term (^ip) is given by 

<f(p) = (A„p“/ 2 ) + ( 7 " - 1 ) P“-" + 


( 2 . 20 ) 

which simplihes to 

Note that 


4^ 

<^{p) = -Ag^{p-^)p-^ + 2p-^. 

d 


dr 

dr‘^ 


P = -r, 


p= -1. 


Observe also that 


and 




9 -1 -2 

jrP = rp , 
or 

^P~' = P~" + 


so that Ag^p =4 + 


. , ^ and 

p 


S{p) = -p Ag^p ) + 2p = 


d-4 

p3 


Observe that in this case, it is reasonable to expect that ijj will be a radial function. 
Assuming so, we can compute the term 

-<^3o(P"^W) = -2rp-V + 


= — 2 rp — p ^ [ 4 ' + 


-2 


d — 1 , / 
— 


= -p M p - ^) V'' - p M + 


= -p 'tp - 

and equation (I2.19p is equivalent to 


d — 5 


d — 1 / 

-V' 


^ -p V' - -p V' 


( 2 . 21 ) 


p v 


-—-p — -p — (d — 4)p ^'ip = 0. 
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We can now verify that ip = 1 + is a solution of the ODE (I2.2ip with the desired 
boundary condition ip\sB = 1- This hnishes the proof of (i). 

We now prove (ii). When d = 4, one can either follow a procedure similar as (i) to 
compute explicitly the solution of the PDE —Ag^r = d — 1 = 3 on 5^. Or denote 
r = |x|, we hnd that a radial solution of —Ag^r = d — 1 = 3 satishes the equation 

6 r 

- 'r(r] —- 'r(r\ = - 

^2 ’ 


- rfr) 


d 2 5 3 


and a solution to this equation satisfying the boundary condition e'^’^gnls^ = ds® is 
given by 


r(x) = log 


1 — \x\ 


1 — \x\ 


Our adapted metric will then be g* = e^'^p~‘^go = g^. We remark alternatively, 

we can hnd the metric g* in dimension d = 4 by a “dimension continuity” argument 
of Branson, by computing the limit 




= lim('0d) 


□ 


Remark 2.4. As pointed out in the introduction, if in (I2.15p we replace the Paneitz 
P 4 operator by the conformal Laplacian and carry out the construction in Proposition 
12.21 ■ we obtain g* = g^. More precisely, if we consider the system 

( 2 . 22 ) = 0 , in 

(2.23) Uf = /, on S^-\ 

(2.24) = 0 on 

OUg* 

where Lg* = —Ag* + then 

(i) when d > 3, g* = go, then ip = Ihj solving the system fl2.22p - fl2.24p 

at s = n/2 -|- 1/2 for the Dirichlet data / = 1). 

(ii) when d = 2 , by solving fl 2 . 12 p for n = 1 on H^, then g* = e^’^gn = go- 


3. Energy identity for the model case in dimension d > 4 

We now write the energy identity in Proposition 12.11 for the metric g* in the model 
case {B‘^, gn) for d > 4 using the results of section [231 Recall that in this 
case the adapted metric is given by g* = ip'^go where '0 = 1-1- ^^p and p{x) = 

Observe that from the dehnition of the energy E[g) in Proposition 12 . 1! as the 
functional obtained by dropping the boundary terms when integrating by parts in 
{{P 4 )gu}udVg we have 

(3.1) Ei{g)[u] = [ u { Pi)gudVg + J^i{g)[u] - J^ 2 {g)[u] - J^ 3 {g)[u], 

Jxd 
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where ^2 and J'j, defined for general u by 


( 3 . 2 ) 

( 3 . 3 ) 

( 3 . 4 ) 


d 




M'^-i 


^2{g)[u] = 

^3{9)[u] = 


f^d—l OtI 


d . . , , 




uh{ng, Vu)dag, 


with h = 4 {Ag — {d — 2)Jg). In view of (13.Ih . let us write the energy identity (I2.14p 
in Proposition 12.11 as 

( 3 . 5 ) f U,(Pt)^U,dV^ + y,(g-)\U,]-y^(g-)\Uf] + y,(g-)\U,] 

(3.6) = 2 (f P^fda - (d - 4) / Qsfda, 

Jsd-l Jsd-1 

where Uj satisfies fl2.15l) - fl2.17p . For the following three lemmas, we will consider 
functions u (which includes the case when u = Uf), which satisfies the following two 
conditions 


( 3 . 7 ) 

( 3 . 8 ) 


u\dB'^ — f, 


d 

dUg* 


u = 0. 


Note that on we have Ug^ = Ug*, thus the Neumann boundary condition fl3.8l) 
can also be written in terms of go. We start by noting the following lemma that holds 
for all dimensions d > 4. 


Lemma 3.1. Let d > 4. For the metric g* and u satisfying (13.7p and (13.Sh we have 
^ 2 .{g*)[u] = 0 for any function f G 

Proof. Recall that J^ 3 {g*)[u] = ® fh{ng*,'Vu)da, where h = 4^4^. — (d — 2)Jg»g* 

Jsd-i 

and by virtue of fl2.17p we have g*{ng*, Vu) = 0, so that 

4 

h{ng*,Vu) = ^RiCg*(ng*, Vu). 

On the other hand, if for d > 4 we write g* = e^'^go where of course 

^ f Pix) for d = 4, 

I ^ log (1 + Vp(^)) ^ ^ 

we have by the transformation law for the Ricci tensor the identity 

RiCg* = (d - 2) - ■^^-^(Ag,r)c/o + dr (g) dr - | Vrp^o^ 

and since r is radial we have 


RiCg» = 6dn (g) dn + 00 go, 
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for radial functions 9 and to and then clearly RiCo*(no*, VnfUm = 0. This implies 
^ 3 (^*)[m] = 0 as needed. □ 

We now relate the boundary terms ^i{g*)[u] and to boundary integrals 

computed in terms of the Euclidean metric Qq. To be more precise, recall that the 
operator (^ 4 )^ satishes the following conformal covariance property: if we write g* = 

4 

then 




and therefore 

/ u{P4)g4u)dVg» 

Jb<^ 

f d+4 . V . V 2d 

= Ui/J ' i -^{ Pi ) g ^{' lljU )' llj <^-^ dVg ^ 

Jb^ 

= / ' lljU { Pi ) g ^{^ l ) U ) dVg ^ 

JB'i 

P d d /* 

= f - -^iu^p)Ag^{u^p)da + / {Ag^'il^ufdVg,, 

Jgd-l aUg^ dUg^ JbA 

= ^i{9o)[iJu] - ^2{9o)[i^u] + [ {Ag^'ipufdVgg, 

J B'^ 


and since / {Ag^ipuYdVg^ is the Paneitz energy of ipu with respect to go , we are 

J B^ 

interested in computing the difference 

d^ii9*)[u] - ^2{g*)[u] - {J^i{go)[ilJu] - J^2igoMu]). 

We will do this not for general u but for u = Uf satisfying fl2.15l) - fl2.17l) . We start 
with the following lemma 

Lemma 3.2. Foru satisfying fl2.15j) - (l2.17p we have the identity 


d^i{9*)[u] - J^i{gQ)[fou] = 


d-4 


/ 


Sd-i 


^ ^90 


u + A/ ) da 


(d-Af 


fda. 


Sd-i 


Proof. On S'® ^ we have Ug* = Ug^ and Ug^ can be identihed with the normal direction 
Since in addition ^i{g*)[u] = 0, we start by noting 


■^i{9\u) - ^i{go)[f}u] = -^i{go)[fju] = 


d 


s<^ 


-1 dn 


■{ipu)Ag^{ipu)da 


go 
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and 


d 


dn 


( 


= + +MV(^9o'0)|r=l, 


d-4 
2 

d-A 


-u 


f 


2 
dr‘i 


d'^ , r \ (d — 4)^ 2 

u\r=i + Am ) + d -— - u^, 




u + A/ ) + d 


from where the result follows. 

We also have 

Lemma 3.3. For any u satisfying (13.7p and (13.Sh . the following identity holds 


^2i9*)[u] - d^2{9o)bl^u] = ^ (f) fAfda + 

gd-l 


d-A 


Proof Since ^/\^\ = ^/\%\ and {g*) = fj ^-4^ \/\%\gQ \ we have 

1 


f^-^uda. 

gd-l OTii^ 


90 


Ag,U = 


\Ai* 


di ( \/W\9'^djU 


2(d — 2)_ ±_, . _ ±_ . 

= ^_4 ^ ‘i-gVlp,Vu)g^+'tp <i-*AggU 

and since fj is radial we obtain 

2{d - 2)_d_ ,du _ ±_ 

and since the normal direction Ug^ coincides with the radial direction A have 


d 


dug ^ 


Ag*u = if <^-4 — 


_2_ d /2(d — 2) ,__d_ gdu ,_ ±_ 

^ A I '■ ' nl\ J A nl\ I nl\ J A 


a-v 4-4 v.-«v.^ + W“A„«|, 


Since u satishes a Neumann boundary condition we have 
d 


dug ^ 


'U I rp —— 


2{d — 2) ,&^u 4 / d . . I 

^ - uu'* + 01^”“ 11-*’ 


and writing 


^ d — 1 d _2 T 

^ (fj'Y' ^ fP CJIT 




we obtain 
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and 


d 


' dr'^ 


-/\g*U 2:'ip |j-=i „M|r’=l “1“ “1“ rN ^go'^\r=l 


A 

dr 


~ d 

4) Q^2 ^ I 4” ^qqU I r=l ■ 

Observe on the other hand that 

(3.9) ^goii^u) = uAg^ij + 2(Vm, VV')go + i^Ag^u, 
and 

2(Vm, VV^)go = 

Since is a constant and u satishes a Neumann boundary condition we obtain 

(3.10) dn^^ {Ag^iiju)) u = ^ -^u + ijAggU^ |^=i, 

,ld‘^U A I ^ A 

(3.11j "ilp Q^2 1^—1 ^u\r=l “h ~Q^^gQU\r=l 


(3.12) 

We have obtained 


3 d — 4 


A I A I 

Au\r=l Ag^u\r=\- 


d 


■{Ag*u)da — 


gd-i dug* ‘9n, 


d 


■{Ag^'ipu)da 


90 


Sd-i 


Sd-I 


di^ ~ d 

u ( -{d - 4)^M|r=i + 2AM|r=i + —Ag^u\r=i ) da 


dr'i 

d'^u 


u ( ~{d - 4)^|.=i - '^^Au\r=i + ^Aggu\r=i ) da 


(f fAfda H- 

2Jsd-i 2 

d 
'2 


2' dr"^ 

5 c /-1 or^ 


a" I . 

M7—Mh=ia(j 


d - A 

fAfda + 

Sd-l ^ 


<92 

fjr^uda. 
gd-i dnt^ 


90 


□ 


We have shown 

Corollary 3.4. For u satisfying fl3.7p and fl3.8p . the following identity holds 

d^ii9*)[u] - J^2{g*)[u] - {J^i{go)[ipu] - J^2{go)[ipu]) 

d(d-4)2 


= -2 


fAfda - 


Sd-i 


rda. 


Sd-i 
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Corollary 3.5. The energy identity fl2.14l) is equivalent to 


(3.13) ^i{go)[^Uf] - = 2 f fPsfda -{d-4) 


Q,fda 


Sd-i 


( 3 . 14 ) 

where Qs is given by ( 12 . 6 p . 


Sd-i 


\Vf\^d<T + 


d{d - Af 


fda, 


Sd-i 


4. Proof of Theorem 

We first assert that there is a class of pseudo differential operators of fractional 
order 27 < d — 1 that are intrinsically dehned on with the conformal covariant 
property. The formula was explicitly computed by Branson |Bra95] . In the special 
case when 27 = 3, we denote the operator as P 3 , the formula of V 3 is given as: 

(4.1) V 3 = {B-1)B{B + 1), 

where B = \J Note that B (and therefore P 3 ) is completely determined 
on spherical harmonics by 

Be = ( 4 + e, 

where is a spherical harmonic of order k. Observe also that for the rough laplacian 
A with respect to the round metric gsd-~^ we have 

Ae^ = -k{k + d-2)i^. 

The operator plays a role in the formulation of the following sharp Moser-Trudinger 
and Sobolev inequalities proved in |Bec93] . Before stating the results in |Bec93] . we 
adopt a convention about normalized measures intended to simplify the proof of 
Lemma [4.41 and Theorem [A1 


Convention 4.1. We will use the following convention for normalized measures: 

• On d(j will be the unnormalized Lebesgue measure with respect to the 

spherical metric gs<i-^- We will use d^ to denote the normalized measure 


( 4 . 2 ) 


d(, 




da 

T(d/2) 


where is the volume of with respect to the metric gs^-i- 

• dx will denote the Lebesgue measure on with respect to the Euclidean 
metric and d/r will be the measure 

(4.3) 

where is the volume of the d-ball B^ with respect to the Euclidean metric, 
note that d|i?'^| = 
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We observe that from the definition of the normalized measures and d/x in fl4.2p 
and (14. 3 p respectively, the measure \S^~^\~^dVgf^ = is such that 




and also 
(4.4) 




\B' 


\S^ 


— = 


d 


We now state the main inequalities in |Bec93] that we are interested in. 

Theorem 4.2 f |Bec93] h Let f G and let V 3 be the operator defined on 

S‘^-^ by dH]). Then 

(a) If d = 4: and (p G C°°{S^) then 


log 




(b) //d > 4 and f G C°°{S'^ then 


ad 


ISd-l 

, d(d-2)(d-A) 

where = ——&—-■ 


2(d-l) 


2{d-l) \ (d-4) 

'^df] < (t mfdf, 

igd-i 


The operator V 3 above stated is closely related to the scattering operator P 3 men¬ 
tioned in Proposition 12.1[ in fact, Branson proved the following 

Theorem 4.3 f [Bra93b] ). On the model space go) for d > 4 and for any g 

in the same conformal class of go, we have 

( 1 ) P^ = Vo. 

( 2 ) 03(g) = PiV:,{l) = 

We will now start the proof of Theorem A by hrst justifying the boundary condition 
fll.Tp imposed in the statement of the theorem. 


Lemma 4.4. A function u satifies (13.7p and 
is in ^f, where denotes the class of functions 


C^^iB'^) :v = f on and 


if and only if the function v = fiu 


dv 


d-4 


dn 


f on S' 


d-l 


go 


Proof of Lemma 4-4 This follows directly from the properties of the function 
which equals one and whose Neumann derivative equals — on the boundary 
of the Euclidean ball B'^. □ 
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Proof of Theorem . HI We start by noting that since Qsig*) = we can rewrite 

the energy identity fl3.13p as 


fPsfda =^J^i{go)['i/jUf] - \^2{go)[ijUf] 

gd—1 Z Z 


Sd-i 


\Vf\^da 


d{d — 4) 


Sd-i 


\f\^da, 


and by Theorem 14.31 we see that 


( 4 . 5 ) 


fP3fd.y = f fVzfda^ 

gd-l 


On the other hand, if we set Vf = iIjU j we see that Vf is a biharmonic extension of / 
to and Vf is in ^f, Since Vf is bi-harmonic, we have 

^d^iigoMUf] - ^J^2{goMUf] = ^J^iigo)[vf] - ^^2(^o)N 


1 

2 


I Bi 


{Ag^Vffdx. 


from Theorem 14.21 part (b), (14. 5 h and 




2(d-l) 


Sd-1 


Ifl^df 


d-4 

2(d-l) 


< 


Sd- 

^.1 


(4.6) 


+ 


we obtain 

fPsfd^ 

{Ag^Vffd^i + 

fdf. 




Sd-i 


|v/pde 


did - 4) 


5 d-i 


I 


{Ag^v^dg. is minimized by 


Thus for a general function n G it is clear that 

bi-harmonic functions in which yields the desired inequality. The assertion on 
the functions for which equality in 11.81 is attained follows from an observation made 
in |Bec93] where it is proved that equality in Theorem 14.21 part (b) is attained for 
functions of the form fiff) = c|l — where c is a constant. We remark that 

for V G ^/, inequality fl4.6p is the same as fll.81) in the statement of Theorem [Al after 
changing the notation dg back to da, dx. □ 


5. Proof of Theorem I 


i-hl= 


Recall that when d = 4, the metric g* is given by g* = e'^^go where p{x) = 
and therefore, for the Paneitz operator P 4 we have 

(6.1) (P4),. = e-‘'’(A„)l 

Fixing a function ip G recall that we have dehned the class given by 

Cg, = Iw e : w\s3 = p>, = 0 
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Assuming that w^p is a bi-harmonic function in we hrst recall that the energy 
identity fl2.14p in the case d = 4 is of the form 

(5.2) Ei{g*)[w^] = 2 I ipPs^pda. 

Js3 

Using the identities fl3.5p and fl3.6p together with the proof of Lemma 13.11 we reach 
the conclusion that fl5.2p is equivalent to 

(5.3) -J^ 2 {g*)[wf] = 2 (f ipPspida. 

Js^ 

We now compute the difference — d^ 2 igo)[wip\- We start by noting that for 

any w in we have 

d d 

d 

= u—(AgoI S3 + 2 Ay?, 

OrigQ 

thus fl5.3p takes the form 


(5.4) 


-^ 2 ( 1 / 0 ) M 


- 2 


® ipAipda = 2 
Js3 


(p ipP:i<pda. 
/S3 


Integrating by parts we obtain 


(5.5) 



(AgoMIg,) dVgg 





We now prove Theorem iBl 



ipP^ipda. 


Proof of Theorem [3 Normalizing the measure da to obtain df and assuming that 
<^ = 0 we see from Theorem 14.21 part (a) that 

log (£= _U_£ |V4=|U5. 

and clearly Wg, minimizes J^ 4 {AggWg,Ydg. in the class and this yields the desired 
inequality. As in the proof of Theorem the assertion about the functions for which 
equality is attained in fll.lOp follows from the form of the functions for which equality 
is achieved in Theorem 14.21 part (a). As proved by Beckner in |Bec93] . all of these 
functions are of the form y? = — log |1 — (zq, 01 + ^ where c is a constant, f ^ and 
Ool < 1 is hxed. □ 


6. Proof of Theorem I 


We start by describing the term l 2 [(?,tc] appearing in fll.l3p explicitly, however, 
before writing out the formula proved in CQ97a| , we need to describe important 
boundary analogues of the Paneitz operator and of Q-curvature for compact manifolds 
with boundary in dimension 4. Consider a manifold with boundary (M^, dM, g) and 
let g be the restriction of g to dM. Let n be the unit outward normal and let us use 
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greek letters to denote directions tangent to the boundary, i.e., Qap 
that with this convention the second fundamental form is given by 


( 6 . 1 ) 


1 d 


Qap- Observe 


and let us use H to denote the mean curvature, i.e., H 
operator dehned in |CQ97a| is 


The boundary 


( 6 . 2 ) 

1 d d 2 1 d 

p\w = —- A-^w - -HAw + + -VH -Vw + iP- 2J)-^w, 

2 on on 3 3 on 


where A, V are the boundary counterparts of Ag, V. Note that this time the function 

—4 

w is dehned in all of B . The convention that we use for the Laplacian is that AgW 
stands for 


AgW = VnVntC + g°‘^Va'^l3W, 

where Vn'^nW means two covariant derivatives with respect to the metric g in the 
normal direction n. The terms F and J and (I6.2jl are used to denote 

• W — IP 

R 

• J = he., the restriction to the boundary of the ambient scalar curvature 

Rg. 

We also have a boundary curvature T 3 given by 


(6.3) t,^-~J + JH- (G, II) + - itr,:(lC + ^AH 

where 


_ pa 


• (G, II) Ran0n^^a0y 
For the model case (5'^, S^, 


we easily see that 


(6.4) II = gs3, 

(6.5) H = 3, 


and {G, II) 
go = ds"^ 

( 6 . 6 ) 

and 

(6.7) 


F = J = 0 so that if n is the outer normal to S'^ with respect to 


1 d 


d 


P^ = ~^Ag,-A^-A 


2 dn 


dn 


Tg = - -tr(Il3) = 2. 

9 3 


The formula proved by Chang and Qing in |CQ97a| is the following 
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( 6 . 8 ) l2[g,w]=h2[go,w] +^D[go,w], 

where 

(6.9) 62 ( 5 ',^ = 7 / w{Pi)gWdVg + \ I wQgdVg + l- I wP^wda+ [ wT^da, 

^ Jb‘^ ^ Jb‘^ ^ Js^ Js^ 

and 0 

(6.10) D[g,w] = — (f {—Rg — + 3F)-^wda — (f H^-^wda — (f HAwda. 

Jg3 3 on Jg3 on Jg3 

We have used in fl6.9p the same notation employed in |CQ97a| . For the model case 
S^,go), a simple computation shows that 

( 6 . 11 ) D[9„w]=3£^(Au,-^u,y 

It is well known that the Paneitz operator {P4)g and Q-curvature {Q4)g have con¬ 
formal covariance properties, but it turns out that Pg^ and Tg are also conformally 
covariant. More precisely, letting gr = e‘^'^g,'we have the transformation laws 

( 1 ) (7)„ = 

( 2 ) (PS),T + (Ts), = (Ts),y\ 

Applying the conformal covariance properties above to the pair of metrics go and 
g* = e'^^go, we easily obtain the following 

Lemma 6.1. The metric g* has the following properties 

(1) (Qi)^* = (^ 4)50 = 

( 2 ) (Ts),. = in),, = 2 , 

hence 

( 6 . 12 ) b2[g*,w] = b2[go,w], 

for all we C°°(P^). 

It was pointed out by Chang and Qing in |CQ97b| that on the model space 
{B^,S^,\dxy) where again go = \dxy is the Euclidean metric, there is a close re¬ 
lationship between the operator P| and Beckner’s Pg operator. In fact, we have 

Lemma 6.2 ( |CQ97b| , Lemma 3.3 and Corollary 3.1). Suppose that w satisfies 

(6.13) Al^w = 0 on 

(6.14) w\s3 = 

Then 

Pgwiss = Vsip. 


^Note that we are using a different sign convention than the one used by Chang and Qing in 
|CQ97a[ CQ97b| since they consider the inner unit normal and not the outer normal. 
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Remark 6.3. We point out that Lemma 16.21 holds regardless of the value of the 
normal derivative of w ai S^. 

As before we will be interested in the class of functions w defined on satisfying 
-^w\s3 = 0. In addition to conformal invariance, the functional I 2 
satishes the following property for the metrics go and g* = e^^go'- 

Lemma 6.4. For any function w in C^p and g* = e^^go where again p = ^ we 

have 


= l2[go,w + p\. 

Proof. Since {B'^,g*) has a totally geodesic boundary and w G we observe that 
D[g*,w] = 0. It then follows from Lemma [6.11 we have 

(6.15) h[9%w] = b2[g\w] = b2[go,w]. 

From the expression for D[go,w] in fl6.1ip we have 


(6.16) 


l2[5'o, w + p]= b 2 [go, w + p] + Y^D[go, w] 

= b 2 [go, w + p] + ^ (f {w + p)nda 

4 Js3 


53 


{w + p) under. 


We now compute all terms involved in the expression of b 2 [go, w + p]. The function p 
satisfies the following equations with respect to the flat metric 

= 0 . 

P3V = 0, 

and from these identities we observe that 


(6.17) 

(6.18) 


'54 


{W+p){PA)go{W + p)dVg^= W{PA)ggWdVg^+ p{PA)goWdVg^, 


I 


IB'^ 


W + p){P3)9oiw + = f w{Ps)g^wda. 

S3 Js^ 


Our goal now is to compute the term p{PA)g^wdVgQ. Integrating by parts we have 


I 


p{PA)ggWdVg^ = I p{Al^)wdV^ 


I B'T 


90 


I B‘r 


i^gop) i^9oW) dVg^ + f Pit i.^9oW) da - 


d 


S3 


dn 


S3 


dp 


= -4 


/B4 


Ag^wdVgo + i AggWda, 


S3 




(6.19) 







SOBOLEV-TRACE INEQUALITIES OF ORDER FOUR 


23 


Using (I6.19P together with the expression + Aw we conclude that 


( 6 . 20 ) 


IB* 


W (P4)g(, wdVgo = f {Wnn “ Wn)da. 


S3 


Combining fl6.17p . fl6.18p . Lemma ICTl and fl6.20p we obtain 


( 6 . 21 ) 


Hdo, W + p]= fealfi'o, w] + ^ f Wnn)da, 

4 JS3 


and combining fl6.2ip with flb.lSp and fl6.16p we conclude that 

h[9*,w] = l2[go,w + p], 


for all w G as needed. 


□ 


We now prove Theorem O 

Proof of Theorem o/O We start by proving part (1). Since w & C^, we have from 
equation fl6.12l) 

h[9*,w] = b2[9o,w]. 

Recall that for the model space S^, go) and for w E we have 

(6.22) {P4)goW = Aj^w, 

Id 

iP3)goW = 


(6.23) 

(6.24) (Tg),, = 2, 

(6.25) (Q4),o = 0, 
and a simple computation shows that 

1 


h[9*,ui] = b2[9o,w] = 


IB* 


{Agf^w)‘^dVgo + - j) \'V(p\'^da+ 2 (j) ipda. 


From Theorem [B] we have the lower bound 


M9*,w] = 


> 


4 , 
471^ 


I (AgfjW) dVgg + 
B* ^ 


S3 


\Vppda + 2 ® If da 
S3 is3 


log '^">da^ f 


and since 


® e^da = |S'^| = 27r^ we obtain 

Js^ 

12 ( 5 '*, w] > + 2 ® fda = —2 (h fda + 2 (h fda = 0 . 


Now, if 12 ( 5 '*, w] = 0 we have 

\L 


1 f ~ 

{Ag^w)^dVg, + - \Vf\^da = — log 


27r2 


S3 


e^^^-Pda^ , 
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and from Theorem [B] it follows that ip{^) = — log |1 — {zq, 01 + c where Zq G is 
hxed, ^ e c is a constant and ta is a bi-harmonic extension of (p in C^. 

Part (2) follows from combining part (1) with Lemma [6.41 □ 
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